Development of a novel mathematical model that explains SARS-CoV-2 infection dynamics in Caco-2 cells.

Mathematical modeling is widely used to study within-host viral dynamics. However, to the best of our knowledge, for the case of SARS-CoV-2 such analyses were mainly conducted with the use of viral load data and for the wild type (WT) variant of the virus. In addition, only few studies analyzed models for in vitro data, which are less noisy and more reproducible. In this work we collected multiple data types for SARS-CoV-2-infected Caco-2 cell lines, including infectious virus titers, measurements of intracellular viral RNA, cell viability data and percentage of infected cells for the WT and Delta variants. We showed that standard models cannot explain some key observations given the absence of cytopathic effect in human cell lines. We propose a novel mathematical model for in vitro SARS-CoV-2 dynamics, which included explicit modeling of intracellular events such as exhaustion of cellular resources required for virus production. The model also explicitly considers innate immune response. The proposed model accurately explained experimental data. Attenuated replication of the Delta variant in Caco-2 cells could be explained by our model on the basis of just two parameters: decreased cell entry rate and increased cytokine production rate.

HLA-A*01:01 allele diminishing in COVID-19 patients population associated with non-structural epitope abundance in CD8+ T-cell repertoire.

In mid-2021, the SARS-CoV-2 Delta variant caused the third wave of the COVID-19 pandemic in several countries worldwide. The pivotal studies were aimed at studying changes in the efficiency of neutralizing antibodies to the spike protein. However, much less attention was paid to the T-cell response and the presentation of virus peptides by MHC-I molecules. In this study, we compared the features of the HLA-I genotype in symptomatic patients with COVID-19 in the first and third waves of the pandemic. As a result, we could identify the diminishing of carriers of the HLA-A*01:01 allele in the third wave and demonstrate the unique properties of this allele. Thus, HLA-A*01:01-binding immunoprevalent epitopes are mostly derived from ORF1ab. A set of epitopes from ORF1ab was tested, and their high immunogenicity was confirmed. Moreover, analysis of the results of single-cell phenotyping of T-cells in recovered patients showed that the predominant phenotype in HLA-A*01:01 carriers is central memory T-cells. The predominance of T-lymphocytes of this phenotype may contribute to forming long-term T-cell immunity in carriers of this allele. Our results can be the basis for highly effective vaccines based on ORF1ab peptides.

On the peptide binding affinity changes in population-specific HLA repertoires to the SARS-CoV-2 variants Delta and Omicron.

OBJECTIVE: To investigate the changes of Spike protein-HLA binding affinity profiles between the Wuhan strain and two dominant variants, the Delta and the Omicron strains, among the Taiwanese, the British and the Russian populations. METHODS: The HLA frequencies and the HLA-peptide binding affinity profiles in the T-CoV database were combined to conduct the study. We focused on the public alleles in the three populations (HLA-A, HLA-B, HLA-C, HLA-DRB1, and/or HLA-DPA1/DPB1 alleles) and the altered peptides of the spike protein (compared to the Wuhan strain) in the Delta G/478K·V1 (B.1.617.2 + AY.1 + AY.2) and the Omicron (BA.1) strains. RESULTS: For the Delta strain, tight bindings of the altered peptides to the HLA alleles decrease in all three populations and almost vanish in the Taiwanese population. For the Omicron strain, tight bindings are mostly preserved for both HLA classes and in the Taiwanese and the British populations, with a slight reduction in HLA class II in the Taiwanese (1.4%), while the Russian population preserves a relatively high fraction of tight bindings for both HLA classes. CONCLUSION: We comprehensively reported the changes in the HLA-associated SARS-CoV-2 Spike protein peptide binding profiles among the Taiwanese, the British, and the Russian populations. Further studies are needed to understand the immunological mechanisms and the clinical value of our findings.

Alterations in SARS-CoV-2 Omicron and Delta peptides presentation by HLA molecules.

The T-cell immune response is a major determinant of effective SARS-CoV-2 clearance. Here, using the recently developed T-CoV bioinformatics pipeline (https://t-cov.hse.ru) we analyzed the peculiarities of the viral peptide presentation for the Omicron, Delta and Wuhan variants of SARS-CoV-2. First, we showed the absence of significant differences in the presentation of SARS-CoV-2-derived peptides by the most frequent HLA class I/II alleles and the corresponding HLA haplotypes. Then, the analysis was limited to the set of peptides originating from the Spike proteins of the considered SARS-CoV-2 variants. The major finding was the destructive effect of the Omicron mutations on PINLVRDLPQGFSAL peptide, which was the only tight binder from the Spike protein for HLA-DRB1*03:01 allele and some associated haplotypes. Specifically, we predicted a dramatical decline in binding affinity of HLA-DRB1*03:01 and this peptide both because of the Omicron BA.1 mutations (N211 deletion, L212I substitution and EPE 212-214 insertion) and the Omicron BA.2 mutations (V213G substitution). The computational prediction was experimentally validated by ELISA with the use of corresponding thioredoxin-fused peptides and recombinant HLA-DR molecules. Another finding was the significant reduction in the number of tightly binding Spike peptides for HLA-B*07:02 HLA class I allele (both for Omicron and Delta variants). Overall, the majority of HLA alleles and haplotypes was not significantly affected by the mutations, suggesting the maintenance of effective T-cell immunity against the Omicron and Delta variants. Finally, we introduced the Omicron variant to T-CoV portal and added the functionality of haplotype-level analysis to it.

Immunogenic epitope panel for accurate detection of non-cross-reactive T cell response to SARS-CoV-2.

The ongoing COVID-19 pandemic calls for more effective diagnostic tools. T cell response assessment serves as an independent indicator of prior COVID-19 exposure while also contributing to a more comprehensive characterization of SARS-CoV-2 immunity. In this study, we systematically assessed the immunogenicity of 118 epitopes with immune cells collected from multiple cohorts of vaccinated, convalescent, healthy unexposed, and SARS-CoV-2-exposed donors. We identified 75 immunogenic epitopes, 24 of which were immunodominant. We further confirmed HLA restriction for 49 epitopes and described association with more than 1 HLA allele for 14 of these. Exclusion of 2 cross-reactive epitopes that generated a response in prepandemic samples left us with a 73-epitope set that offered excellent diagnostic specificity without losing sensitivity compared with full-length antigens, and this evoked a robust cross-reactive response. We subsequently incorporated this set of epitopes into an in vitro diagnostic Corona-T-test, which achieved a diagnostic accuracy of 95% in a clinical trial. In a cohort of asymptomatic seronegative individuals with a history of prolonged SARS-CoV-2 exposure, we observed a complete absence of T cell response to our epitope panel. In combination with strong reactivity to full-length antigens, this suggests that a cross-reactive response might protect these individuals.

Endocytosis and Transcytosis of SARS-CoV-2 Across the Intestinal Epithelium and Other Tissue Barriers.

Since 2003, the world has been confronted with three new betacoronaviruses that cause human respiratory infections: SARS-CoV, which causes severe acute respiratory syndrome (SARS), MERS-CoV, which causes Middle East respiratory syndrome (MERS), and SARS-CoV-2, which causes Coronavirus Disease 2019 (COVID-19). The mechanisms of coronavirus transmission and dissemination in the human body determine the diagnostic and therapeutic strategies. An important problem is the possibility that viral particles overcome tissue barriers such as the intestine, respiratory tract, blood-brain barrier, and placenta. In this work, we will 1) consider the issue of endocytosis and the possibility of transcytosis and paracellular trafficking of coronaviruses across tissue barriers with an emphasis on the intestinal epithelium; 2) discuss the possibility of antibody-mediated transcytosis of opsonized viruses due to complexes of immunoglobulins with their receptors; 3) assess the possibility of the virus transfer into extracellular vesicles during intracellular transport; and 4) describe the clinical significance of these processes. Models of the intestinal epithelium and other barrier tissues for in vitro transcytosis studies will also be briefly characterized.

T-CoV: a comprehensive portal of HLA-peptide interactions affected by SARS-CoV-2 mutations.

Rapidly appearing SARS-CoV-2 mutations can affect T cell epitopes, which can help the virus to evade either CD8 or CD4 T-cell responses. We developed T-cell COVID-19 Atlas (T-CoV, https://t-cov.hse.ru) - the comprehensive web portal, which allows one to analyze how SARS-CoV-2 mutations alter the presentation of viral peptides by HLA molecules. The data are presented for common virus variants and the most frequent HLA class I and class II alleles. Binding affinities of HLA molecules and viral peptides were assessed with accurate in silico methods. The obtained results highlight the importance of taking HLA alleles diversity into account: mutation-mediated alterations in HLA-peptide interactions were highly dependent on HLA alleles. For example, we found that the essential number of peptides tightly bound to HLA-B*07:02 in the reference Wuhan variant ceased to be tight binders for the Indian (Delta) and the UK (Alpha) variants. In summary, we believe that T-CoV will help researchers and clinicians to predict the susceptibility of individuals with different HLA genotypes to infection with variants of SARS-CoV-2 and/or forecast its severity.

Association of HLA Class I Genotypes With Severity of Coronavirus Disease-19.

Human leukocyte antigen (HLA) class I molecules play a crucial role in the development of a specific immune response to viral infections by presenting viral peptides at the cell surface where they will be further recognized by T cells. In the present manuscript, we explored whether HLA class I genotypes can be associated with the critical course of Coronavirus Disease-19 by searching possible connections between genotypes of deceased patients and their age at death. HLA-A, HLA-B, and HLA-C genotypes of n = 111 deceased patients with COVID-19 (Moscow, Russia) and n = 428 volunteers were identified with next-generation sequencing. Deceased patients were split into two groups according to age at the time of death: n = 26 adult patients aged below 60 and n = 85 elderly patients over 60. With the use of HLA class I genotypes, we developed a risk score (RS) which was associated with the ability to present severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) peptides by the HLA class I molecule set of an individual. The resulting RS was significantly higher in the group of deceased adults compared to elderly adults [p = 0.00348, area under the receiver operating characteristic curve (AUC ROC = 0.68)]. In particular, presence of HLA-A*01:01 allele was associated with high risk, while HLA-A*02:01 and HLA-A*03:01 mainly contributed to low risk. The analysis of patients with homozygosity strongly highlighted these results: homozygosity by HLA-A*01:01 accompanied early deaths, while only one HLA-A*02:01 homozygote died before 60 years of age. Application of the constructed RS model to an independent Spanish patients cohort (n = 45) revealed that the score was also associated with the severity of the disease. The obtained results suggest the important role of HLA class I peptide presentation in the development of a specific immune response to COVID-19.

HIF Prolyl Hydroxylase Inhibitors for COVID-19 Treatment: Pros and Cons.

The review analyzes the potential advantages and problems associated with using HIF prolyl hydroxylase inhibitors as a treatment for COVID-19. HIF prolyl hydroxylase inhibitors are known to boost endogenous erythropoietin (Epo) and activate erythropoiesis by stabilizing and activating the hypoxia inducible factor (HIF). Recombinant Epo treatment has anti-inflammatory and healing properties, and thus, very likely, will be beneficial for moderate to severe cases of COVID-19. However, HIF PHD inhibition may have a significantly broader effect, in addition to stimulating the endogenous Epo production. The analysis of HIF target genes reveals that some HIF-targets, such as furin, could play a negative role with respect to viral entry. On the other hand, HIF prolyl hydroxylase inhibitors counteract ferroptosis, the process recently implicated in vessel damage during the later stages of COVID-19. Therefore, HIF prolyl hydroxylase inhibitors may serve as a promising treatment of COVID-19 complications, but they are unlikely to aid in the prevention of the initial stages of infection.

Potential role of cellular miRNAs in coronavirus-host interplay.

Host miRNAs are known as important regulators of virus replication and pathogenesis. They can interact with various viruses through several possible mechanisms including direct binding of viral RNA. Identification of human miRNAs involved in coronavirus-host interplay becomes important due to the ongoing COVID-19 pandemic. In this article we performed computational prediction of high-confidence direct interactions between miRNAs and seven human coronavirus RNAs. As a result, we identified six miRNAs (miR-21-3p, miR-195-5p, miR-16-5p, miR-3065-5p, miR-424-5p and miR-421) with high binding probability across all analyzed viruses. Further bioinformatic analysis of binding sites revealed high conservativity of miRNA binding regions within RNAs of human coronaviruses and their strains. In order to discover the entire miRNA-virus interplay we further analyzed lungs miRNome of SARS-CoV infected mice using publicly available miRNA sequencing data. We found that miRNA miR-21-3p has the largest probability of binding the human coronavirus RNAs and being dramatically up-regulated in mouse lungs during infection induced by SARS-CoV.

Integrative analysis of miRNA and mRNA sequencing data reveals potential regulatory mechanisms of ACE2 and TMPRSS2.

Development of novel approaches for regulating the expression of angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) is becoming increasingly important within the context of the ongoing COVID-19 pandemic since these enzymes play a crucial role in cell infection. In this work we searched for putative ACE2 and TMPRSS2 expression regulation networks mediated by various miRNA isoforms (isomiR) across different human organs using publicly available paired miRNA/mRNA-sequencing data from The Cancer Genome Atlas (TCGA) project. As a result, we identified several miRNA families targeting ACE2 and TMPRSS2 genes in multiple tissues. In particular, we found that lysine-specific demethylase 5B (JARID1B), encoded by the KDM5B gene, can indirectly affect ACE2 / TMPRSS2 expression by repressing transcription of hsa-let-7e / hsa-mir-125a and hsa-mir-141 / hsa-miR-200 miRNA families which are targeting these genes.